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Abstract. A number of different methods have been proposed to identify 
' unanticipated burst sources of gravitational waves in data arising from LIGO and 

other gravitational wave detectors. When confronted with such a wide variety of 
methods one is moved to ask if they are all necessary: i.e., given detector data that 
^ ' is assumed to have no gravitational wave signals present, do they generally identify 

, the same events with the same efficiency, or do they each 'see' different things in the 

detector? Here we consider three different methods, which have been used within the 
. LIGO Scientific Collaboration as part of its search for unanticipated gravitational wave 

I bursts. We find that each of these three different methods developed for identifying 

candidate gravitational wave burst sources are, in fact, attuned to significantly different 
features in detector data, suggesting that they may provide largely independent lists 
of candidate gravitational wave burst events. 



o 



PACS numbers: 04.80.Nn, 95.55.Ym, 07.05.Kf 



X, 

^ , 1. Introduction 



The great adventure of gravitational wave detection is the search for the unknown: 
bursts of gravitational waves from sources never anticipated. Searching for the known is 
a different kind of endeavor than searching for the unknown. In the former case, one is 
asking an affirmative question of the data-e.g., is a signal like the one I know arises from 
an inspiraling compact binary system present in the data? In the latter case, however, 
the best one can do in identifying candidate unanticipated bursts is ask if the data now 
is significantly different than the long-run character of the noise. 

The gravitational wave detection community has developed a number of different 
methods whose goal is to look for such anomalous behavior in the gravitational wave 
detector data [H 12 El IH Ej • When confronted with such a wide variety of methods 
one is moved to ask if they are all necessary: i.e., do they generally identify the same 
events with the same efficiency, or do they each 'see' different things in the detector? 
The LIGO Scientific Collaboration has used several in its own internal analyses but 
few burst gravitational wave data analysis method comparisons have been performed. 
Two prominent comparisons have been published [3 |H] but neither of these investigated 
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the methods on actual detector data or focused their observations on the results of 
each method using the current practical application of each method. Here we consider 
three different analysis methods-SLOPE j3], BlockNormal jS], and Q Pipeline [Hlini^and 
examine whether they are more or less sensitive to the same types of burst events in the 
detector data using the current practical tunings of each ETC. 

In Section|21we describe the three different analysis methods-SLOPE, BlockNormal 
and Q Pipeline-that we evaluate in this work. Our discussion of SLOPE is longer than 
that of BlockNormal or Q Pipeline because references exist for the implementation of 
these, but not for the implementation of SLOPE that we use here. In Section El we 
look in detail at the response of all three to 'pure' noise, comparing the false events 
identified for each and looking for correlations among them. Finally, we summarize our 
conclusions in Section 0] 

2. Several Burst Gravitational Wave Data Analysis Methods 

Every analysis pipeline contains a component that performs the first identification of 
what may later become a candidate gravitational wave event. In the LIGO Scientific 
Collaboration this pipeline component is referred to as an event trigger generator, or 
ETG. While later stages of the pipeline may prune the list of candidate events, it is the 
ETG that determines the initial list of possibilities. We focus our attention on the three 
ETGs referred to as SLOPE, BlockNormal and Q Pipeline and the list of triggers each 
generates. In this section we briefly describe each of these ETGs. 

2.1. SLOPE 

As its name suggests, the SLOPE ETG looks for segments of data that show a mean 
trend: i.e., that show a statistically significant slope that persists over some period of 
time. White noise will show, in the mean, no slope, while a gravitational wave burst 
incident on the detector will show at the very least a rise and fall in amplitude associated 
with the beginning and end of the burst. SLOPE was originally proposed by Pradier, 
et al. |lj and was used internally by LIGO as part of its initial data analysis efforts |lUj . 

SLOPE has previously been applied to LIGO data during the first LIGO science 
data run in 2002. Ultimately, unfortunate choices in SLOPE'S application, as well as 
immature data conditioning, led to SLOPE being excluded from the burst gravitational 
wave upper limit calculation ^0]- Described here are the two major modifications to 
SLOPE from its first LIGO application 

2.1.1. Thresholding In the absence of a gravitational wave signal, assuming that the 
detector noise is Gaussian and white, the probability of measuring a slope of magnitude 
greater than \m\ is: 




(1) 
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where /z^ and cr^ (the mean and variance of the slopes produced by SLOPE) are given 
by Equations 2.4 and 2.5 in jjj: 

2 _ l^/dafg 2 

/im = (3) 

and where cr'^ata is the variance of the data being analyzed by SLOPE. These also 
approximate the data as Gaussian. 

The measured value of \m\ is compared against a threshold mthr-esh > 0. nithresh is 
determined based on the probability E using the cumulative distribution function C: 
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where C is the probability of measuring a slope between — oo and nithresh and P is the 
normal probability distribution of slope measurements given a^. and Since the mean 
of this distribution is zero (since = 0), P is an even function with respect to the 
measured slopes and comparing the absolute values of the slopes to the threshold slope 
is justified. S is identified as the confidence that a measured slope is not accidental: 

S = 2C - 1 = erf (5) 

\ 0'mV2 / 

Solving for \mthresh\ yields the expression for the threshold slope: 

\mhTesh\ = CTmV2eTi~^{T.) (6) 

If for some windowed data \m\ > nithresh then a potential gravitational wave candidate 
is identified. The value of nithresh is chosen to ensure that the rate of false candidate 
identifications is an acceptable value. An illustration of this thresholding method is 
shown is Figure [H 



2.1.2. Temporal Slope Clustering When applying successive windows to the data, a 
phenomenon we refer to as mult i- windowing occurs. Multi-windowing occurs when a 
single event spans multiple windows, leading to multiple triggers. Before the SLOPE 
ETC reports a candidate event triggers that are separated by less than a window 
duration are clustered together into a single event, whose duration is reported as 
appropriately greater than a window length. 



2.2. BlockNormal 

BlockNormal is the name given to an ETG that looks for change-points in the statistics 
of the timeseries data: i.e., points where the mean statistics-in the case of BlockNormal, 
the data mean and variance-changes. Gravitational waves are uncorrelated with 
detector noise and, thus, will lead to a change in both mean and variance; so, change- 
points in the noise statistics mark epochs in which evidence for a gravitational wave 
burst may be present in the data. An overview of the entire BlockNormal data analysis 
pipeline is found in jS]. 
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Figure 1. Illustration of SLOPE thresholding on a timeseries (dotted line). Any 
measured slopes (solid gray line) that are greater than the absolute value of the 
threshold slope (above the top black dashed line or below the bottom dashed black 
line) are marked as triggers and later clustered together. 



Once BlockNormal identifies a set of cliange-points in a timeseries, tlie mean 
power in eacli block (i.e., tlie epocli between successive cliange-points) is calculated and 
compared to a threshold. If the power exceeds the threshold then that block, clustered 
together with any adjacent blocks whose mean power also exceeds the threshold, is 
reported ClS db candidate gravitational wave event. 



2.3. Q Pipeline 

The Q Pipeline ETG is a multi-resolution time-frequency search for excess power in 
the detector output. First, the data are whitened using zero phase linear prediction 
and projected onto bases that are logarithmically spaced in frequency and Q (quality 
factor), and linearly in time. Significant tiles whose amplitude significantly exceeds 
that expected for white, Gaussian noise are reported as triggers assuming white noise 
statistics and the most significant set of non-overlapping tiles are reported [HI E] ■ 
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3. Strongest Accidental Event Characterization 

Our principal objective is to determine the degree to which each ETG sees the data 
differently. One straightforward way of addressing this question is to run each over the 
same data set, which may be real detector data or simulated noise, and ask whether the 
events that are reported by each are the same or different. 

3.1. Methodology and Results 

We ran each ETG over the same set of data, drawn from the LIGO S3 science run, 
with the thresholds of each set so that they all identified 'candidate' events (really, 
false alarms) at the same rate and with the same mean duration. Each ETG reported 
back a list of events, each associated with an ETG-dependent amplitude, start-time and 
duration. These three event lists form the basis for our comparison of the three ETGs. 

The event amplitude reported by each ETG depends on the character of the ETG; 
however, for identical events the amplitude increases monotonically with the event 
energy. The events in each list were assigned a rank in their respective list, determined 
by the reported amplitude of each. If the ETGs are similar we expect that the top- 
ranked events in each list should be correlated with each other: i.e., the highest ranked 
events across each list should have similar start times, the second highest ranked events 
also similar start times, and so on. 

Figure El shows the start times and amplitudes of the ten strongest events from each 
ETG from the list generated in one 600 second interval of data. In this segment of data 
it is clear that some events-for example, the event at GPS time 75198177-are identified 
by all ETGs; however, the ranks associated with the events range considerably and, 
while they are closely spaced in time, a more detailed examination of the data shows 
that there is no overlap of events between the ETGs. This suggests that while there 
is likely some noise artifact in this area that caused this cluster of triggers, the ETGs 
did not process it equivalently. Additionally, it is clear that some ETGs identify strong 
events at places where others do not. 

Together, these two observations strongly suggest that each ETG 'sees' the data 
differently. 

3.2. Rank Comparison of Triggers Identified by Multiple ETG's 

Since there were no obvious similarities in the temporal location of the strongest triggers 
between ETGs and no specific obvious similarities or differences in the timeseries 
properties upon inspection of the strongest ETG triggers, we looked to see if there was a 
correlation between the triggers that were detected by ETGs. Focusing events identified 
by each ETG in 600 second duration data segments we identified those candidates that 
were reported by all three ETGs and that shared at least one data sample in common. 
Each of these triply-coincident triggers was then given three ranks, one associated with 
its relative amplitude in each ETG. 
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Figure 2. Timeline of the top 10 strongest triggers from a single data segment for BiockNormai, SLOPE and Q Pipeline. 
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Rank/Rank Plot of All BlockNormal/SLOPE/Q Pipeline Coincident triggers 
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Figure 3. Scatter plot showing the rank of a triply coincident trigger in the 
BlockNormal ETG to the rank as detected in the SLOPE ETC. 

There are 253 triply-coincident triggers which are 2.057% of the triggers admitted 
by each ETG (top 100 triggers from 123 data sets). However, there is a muhi-triggering 
issue with Q Pipehne; that is, since Q Pipehne produces single points where the trigger is 
localized, there is a higher probability that there can be more than one Q Pipeline trigger 
contained in either a BlockNormal or SLOPE trigger. There are 233 unique triply- 
coincident triggers between BlockNormal and SLOPE and 253 unique triply-coincident 
triggers between BlockNormal and Q Pipeline, and SLOPE and Q Pipeline. Regardless, 
the small percentage of triply-coincident triggers compared to the overall number of 
triggers from each ETG reinforces that there is not an equivalence between the ETG in 
the triggers that they all identify. 

We then plotted the ranks against each other in a scatter plot, as in Figure 01 If 
the two ETGs identified events identically, or near identically, we would expect a strong 
correlation among the ranks, which would show up as a clustering in the scatter plot 
along a diagonal line radiating from the origin. By eye. Figure 0] shows evidence for a 
very weak correlation in the events identified by BlockNormal and SLOPE, but this is 
not a correlation in rank equivalence between ETGs. Similarly, Figure |S1 shows evidence 
for a somewhat stronger correlation in the events identified by jointly by SLOPE and 
Q Pipeline, but not in rank equivalence. Figure |3] shows (again, to the eye) no evidence 
for a correlation among events identified jointly by BlockNormal and Q Pipeline. 
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Rank/Rank Plot of All BlockNormal/SLOPE/Q Pipeline Coincident triggers 
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Figure 4. Scatter plot showing the rank of a triply coincident trigger in the 
BiockNormal ETG to the rank as detected in the Q Pipeline ETC 

4. Conclusions 

This investigation sought to determine if several LIGO burst ETGs revealed similar or 
different information about the data even though each method processed the data in 
a different way. Initial inspection of the temporal locations of the strongest triggers 
in individual data sets showed that the events identified by each ETG do not form a 
strongly overlapping set and that, even when we focus just on events that are identified 
by all three ETGs, each ETG assigns a different relative significance to these events. 
We conclude that these different ETGs view the data from different perspectives and 
that candidate gravitational wave events identified in one may be missed by others. 
Correspondingly, we recommend that multiple methods of searching for gravitational 
wave bursts continue to be pursued and used in gravitational wave data analysis efforts. 

Acknowledgments 

Special thanks to Shurov Chatterji for helping to incorporate Q Pipeline into this study 
and to John McNabb for his continuing input on comparison methodology. We are also 
grateful to Shantanu Desai, Keith Thorne, and Peter Saulson for comments and useful 
discussions. 

This work was supported by the Center for Gravitational Wave Physics, and the 



9 



Rank/Rank Plot of All BlockNormal/SLOPE/Q Pipeline Coincident Triggers 

100 I 1 1 1 1 1 1 1 n 

90 - ' . 

80 - _ " . 

70 - ' ■ 

60 - 



O • . . .... 

_i . 

(73 . . . 

40 - ■ ■ 

30- • . • _ ■ . - 

20 - • ■ . • • •■ . ■ _ 

10 . ■ 

Q llI^ : LI '. • I •' '—I l! i ! I • u ; !_i : i !j 

10 20 30 40 50 60 70 80 90 100 

Q Pipeline Rank 

Figure 5. Scatter plot showing the rank of a triply coincident trigger in the SLOPE 
ETG to the rank as detected in the Q Pipeline ETC 

National Science Foundation under award PHY 00-99559. The Center for Gravitational 
Wave Physics is supported by the National Science Foundation under cooperative 
agreement PHY 01- 14375. 

References 

[1] W. G. Anderson, P. R. Brady, J. D. E. Creighton, and E. E. Flanagan, "An Excess power statistic 
for detection of burst sources of gravitational radiation," Phys. Rev. D, vol. 63, no. 042003, 2001. 

[2] S. Klimcnko, I. Yakushin, M. Rakhmanov, and G. Mitselmakher, "Performance of the WaveBurst 
algorithm on LIGO data," vol. 21, 2004, pp. S1685-S1694. 

[3] J. Sylvestre, "Time-frequency detection algorithm for gravitational wave bursts," Phys. Rev. D, 
vol. 66, no. 102004, 2002. 

[4] T. Pradier, et al., "Efficient filter for detecting gravitational wave bursts in interfcrometric 

detectors," Phys. Rev. D, vol. 63, no. 042002, 2001. 
[5] J. W. C. McNabb, M. Ashley, L. S. Finn, E. Rotthoff, A. Stuvcr, T. Summerscales, P. Sutton, 

M. Tibbits, K. Thorne, and K. Zaleski, "Overview of the BlockNormal event trigger generator," 

vol. 21, 2004, pp. S1705-S1710. 
[6] S. K. Chatterji, "The search for gravitational wave bursts in data from the second LIGO science 

run," Ph.D. dissertation, Massachusetts Institute of Technology, 2005. 
[7] L. Blackburn, et al., "A first comparison of search methods for gravitational wave bursts using 

LIGO and Virgo simulated data," vol. 22, 2005, pp. S1293-S1302. 
[8] N. Arnaud, M. Barsuglia, M.-A. Bizouard, V. Brisson, F. Cavalier, M. Davier, S. Hello Patrice, 



10 



Kreckelbergh, E. K. Porter, and T. Pradier, "Comparison of filters for detecting gravitational 

wave bursts in interferometric detectors," Phys. Rev. D, vol. 67, no. 062004, 2003. 
[9] S. Chatterji, L. Blackburn, G. Martin, and E. Katsavounidis, "Multiresolution techniques for the 

detection of gravitational-wave bursts," Class. Quant. Grav., vol. 21, pp. S1809-S1818, 2004. 
[10] B. Abbott, et al., "First upper limits from LIGO on gravitational wave bursts," Phys. Rev. D, 

vol. 69, no. 102001, 2004. 

[11] A. L. Stuver "Burst Gravitational Wave Data Analysis Methods: Design, Development and 
Comparison," PhD dissertation, The Pennsylvania State University, 2006. 



